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Abstract 

We study the spectrum of hght flavor baryons in a quark-model framework 
by taking into account the order O(a^) hyperflne interactions due to two- 
gluon exchange between quarks. The calculated spectrum agree better with 
the experimental data than the results from hyperflne interactions with only 
one-gluon exchange. It is also shown that two-gluon exchange hyperflne in- 
teractions bring a significantly improved correction to the Gell-Mann-Okubo 
mass formula. Two-gluon exchange corrections on baryon excitations (includ- 
ing negative parity baryons) are also briefly discussed. 
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1. Introduction 



A number of approaches have been developed for describing the baryon 
spectrum since the 1960^'s^such as the SU(6) model [ll|2|, the quark model 0, 

the Skyrme model [sl,!^], the large Nc baryon 




li] , the bai 

model [Id, lUl, ll2|, ll3|, [JJ, M, \m and so on. These models incorporate 
partly the dynamics of quantum chromodynamics (QCD) and have arrived 
at remarkable (though not quantitatively perfect) success. Recent lattice cal- 
culations have made si gnifi cant progress on the spectrum of bound states and 

offer a bright promise [B, E, 0, m 

|2l[. However, the quark model still re- 
mains a basic and indispensable tool for understanding hadron spectroscopy 
due to its intuition and simplicity as a guide-line to other approaches. The 
quark model has been proved quite fruitful on the study of baryon spectrum, 
decays and moments. In the quark model with dynamics, the spectrum of 
hadrons is dominated by two ingredients: the long range forces to bind quarks 
as hadrons and the short range forces expected from gluon exchanges. The 
long range forces are the confining forces, which are fiavor and spin indepen- 
dent. The short range forces include a Coulomb term, hyperfine interaction 
including the spin-orbit interaction. Actually spin-spin part of the hyperfine 
interaction is the most important short range force that is responsible for 
such prominent features such as the A-A^ and p-vr mass splittings [3]. 

The hyperfine interaction employed in quark model is usually derived 
from one-gluon exchange and it contains a spin-spin and a tensor part. In 
baryons the spin-orbit interaction can be neglected. The next step for a 
further improvement of the hyperfine interaction is to consider not only one- 
gluon exchange but also higher order terms such as two-gluon exchange. In 
fact, such a potential with two-gluon exchange has already been calculated by 
Gupta and Radford |22l. |23| . Quarkonium spectra were investigated with the 



higher order potential [2J, |25|, |26j and the theoretical results agree excellently 



with experiments. Inspired by these results, we try to investigate the light 
fiavor baryon spectrum with the hyperfine interaction of higher order O(a^). 
Since the nonrelativistic quark model is the simplest and most economical 



quark potential model with considerable phenomenological success [2 



29, M 31, 32, 33 



we try to study the nonrelativistic quark model with the 
hyperfine interaction of higher order O(a^). Furthermore, we particularly 
focus our attention on the ground-state baryon masses, as all wave functions 
of ground-state baryons have zero orbital angular momentum so that we 
can neglect the spin-orbit forces which should be better to be treated in a 



2 



relativistic framework. 

In this paper, we adopt the simple harmonic-oscillator quark model with 
the interactions between quarks being harmonic springs and hyperfine inter- 
actions. The hyperfine interaction includes not only one-gluon exchange of 
order a^, but also two-gluon exchange of order . After choosing the tra- 
ditional zero-order wave functions of the SU(6) multiplets (56,0"^), (70, 
(56', 0+), (70, 0+), (56, 2+), (70, 2+) and (20, 1+) at = 0, = 1 and = 2 
respectively, we employ the Hamiltonian and the baryon wave functions to 
calculate the matrix elements. Since the hyperfine matrix elements of 0{as) 
have been calculated in previous studies, we only calculate the hyperfine ma- 
trix elements of 0(af). We then diagonalize the complete matrices to get the 
baryon masses. We find that the masses are mostly in agreement with the 
experimental masses. In addition, the masses of the ground-state baryons 
are better than those with only one-gluon exchange. We also find that the 
two-gluon exchange interactions bring a significant improved correction to 
the Gell-Mann-Okubo mass formula of the ground-state baryons. We also 
briefiy discuss the two-gluon exchange corrections on negative-parity baryons 
and other baryon excitations. 

2. The Hamiltonian and its solutions 

The Hamiltonian in the quark model 3^ is 



H = + Ho + Hhyp, (1) 



i<j 



i<j 

where Vj^^j is the spin-independent potential and -f^hyp the hyperfine in- 
teraction. The potential V^^^j is 

1 



Cnf = o^< + f^M, (4) 
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where anharmonicity U{rij) is a function which only depends on rjj, and the 
exphcit form of U{rij) is not needed in this paper. 

^hyp 



The hyperfine interaction H^^ is 



(5) 



where we add the order a1 interaction. 

The hyperfine interaction of order as derived from the one-gluon exchange 
process is 



^hyp(as 



2a, ,87r 



Sminij 3 
+ 



Si ■ Sj5^{rij^ 



1 ^ 3Si ■ Tij Sj ■ Tjj ^ 



(6) 



where and Si are the mass and spin of the ith quark and fij is the relative 
position of the pair ij of quarks. The first term is called the Fermi contact 
term and operates when the quark pair has zero orbital angular momentum. 
The second term is called the tensor term and operates in states with nonzero 
orbital angular momentum between the quarks i and j. 



The Fermi contact term and the tensor term of order a^ 22|, |2J, |26 
derived from two-gluon exchange are 



-'-'contact ' 



^^S.-S,{[^{26 + 9ln2)]S'{r-,) 

a. 



247r2 



(33-2n,)^V"(^^^^-^-^+^^] 



+ 



21a. 



tensor > 



167r2 
2a 



\n{^/m~m~rij) + 7e 



(7) 



1 r3{Si ■ rij){Sj ■ Tij) 



Sminij rf 



Si ■ Sj] 



"J ' ij 

a. 



Aa ^ 

+ F^(22 " 2n/)[ln(/XGRry) + 7e - ^] 

OTT OTT 6 



2>a. 



IT 



[H^rmmjr.j) + 7e - -]}, 



(8) 



where /xgr is a renormalization scale, the subscript GR refers to the renor- 
malization scheme 34] and nj is the number of effective quark flavors. 
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We then write approximate solutions 30|, l31| by perturbation theory in 
U and Hhyp. In the U = -ffhyp = hmit, Hq becomes 

Ho^K = ^ + ^ + ^ip' + \'), (9) 



where 



2mp 2m\ 2 



P = -^(ri-rs), pp = mp—, (10) 
A = -^(fi + f2 - 2f3), p^ = mx^, (11) 



with 



rup = m, mx = - ■ r, x = md/ms. (12) 

2m + m' 

In the = sector, m = m' = m^; in the 5* = — 1 sector, m = m^, 
m' = rris] in the S = —2 sector, m = rris, m' = 171^ and in the 5* = —3 sector, 
m = m' = rris. 

The solutions to H'^ are wave functions as follows. The wave function 
with iV = is 

3/2 3/2 

^00 = ^^exp[-iay-ia^A^)]. (13) 
The wave functions with = 1 and N = 2 are 

3/2 5/2 -, , 

= ^^A+exp[-^ay-^a^A^)], (14) 

5/2 3/2 . , 

<i = ^^p+exp[-^«y-i«lA^)], (15) 

/;T 3/2 7/2 o 11 

- ^«-)exp[-i«y - lalX% (16) 

9 5/2 5/2 -, . 

< = ;^^^P-Aexp[-i«y-ia^A^)], (17) 

/;T 7/2 3/2 o 11 

- ^«;')exp[-^«y - lalX% (18) 
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1 3/2 7/2 . . 

V^2^2^ = ^^^A^A^exp[-iay-ialA^)], (19) 

5/2 5/2 . , 

< = ^^P+A+exp[-iay-i«^A^)], (20) 

1 7/2 3/2 . , 

^2^2^ = ;7|^^P+P+exp[-iay-iaiA^)], (21) 

5/2 5/2 , . 

< = ^^-^(P+A3-p3A+)exp[--ay--a^A^)], (22) 



where we have shown only the highest state of an orbital angular momentum 
multiplet, which has energy (2np + lp + 3/2)up + {2nx + + ?>/2)uJx and 



a 



p 



{3KmpY/\ ax = {3Kmx)'/\ (23) 
(3K/m,)i/^ ujx = {3K/mxf'^. (24) 



When U differs from zero, it can be shown that in the S* = sector, the 
confinement energies can be determined by three constants Eq, a, and b, 

E[iIjoo] = Eo, (25) 
E[^|J^m] = Eo + a, (26) 

= i^[^oo] =Eo + 2a- h, (27) 

E[roo] = Eo + 2a-h, (28) 

= i?[<J = Eo + 2a- ^b, (29) 

E[^2l] = Eo + 2a-i6, (30) 
E[riL]=Eo + 2a, (31) 
(Clt^l^oo) = -\b, (32) 

(V^2^^;|f/|V^2^) = (33) 

To break SU(3), p and A excitation energies are decreased by (md/mpY^'^ 
and {jrid/mxY^'^ respectively, as they would be in the harmonic limit. We 
use the confinement energies shown in Table [H obtained by considering U 
effects in Refs. 31, 32 
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We then calculate the hyperfine matrix elements. In the nonstrange sec- 
tor, the states are completely symmetric in flavor, spin and space under 
interchange of any two quarks, so it follows that 

{a\H^,,\P)=3{a\Hl%\P), (34) 

where \a) and are any two states. This trick simplifies the calculation. 
In the S = —1 sector, since the states are always symmetric under exchange 
of quarks one and two, then 

{a\Hl^^\P) = (a|ifgp|/3), (35) 

and the hyperfine matrix elements become 

{a\H^y,\/3) = {a\Hl% + 2Hl^^\f3). (36) 

Calculation of the -f^hyp(^i2) rnatrix elements is straightforward. To perform 
the -f^hyp(^i3) matrix elements, we relate to the interaction between quarks 
one and two by using permutations: 

(«K'pl/5) = (a|(23)(23)i/^Jp(23)(23)|/?), (37) 

with 



(23X^p(23) = -HS„ (38) 
(23)|a) = 5^c„,|«'), (39) 



where Ca' is the coefficient of the \a') states determined by spin and space 
functions permutations. The details of the hyperfine matrix elements calcu- 
lations are discussed in Appendix \M The hyperfine matrix elements of the 
S = —2 and S = —3 sectors can be obtained from the S = and S = —1 
sectors by making the interchange rriu mg everywhere. Since the hyperfine 
matrix elements of 0{as) have been calculated 29, 3^, 31|, we only calculate 



the hyperfine matrix elements of ©(a^) derived from two-gluon exchange. 
The results are listed in Tables [TT] and [T21 The detailed calculation of the 
hyperfine matrix elements of higher order are discussed in Appendix [Bl We 
also express the hyperfine matrix elements with O(a^) here in the unit of 
(5, which is 4as«^/3v^27rm^. In this paper, we assume that this value differs 
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slightly 30|, l3l|], which is 270 MeV for considering the results of order ef- 
fects and mixing between the ground states with = and the exited states 
with N = 2. It shows that the values of matrix elements of higher order are 
in the range (0.0014-0.0334)5, i.e, the values of the matrix elements are in the 
range 0.4-10 MeV, which are one order less than those 2-170 MeV obtained 
by order a^. At last we diagonalize the complete Hamiltonian matrices to 
get the baryon masses. The hyperfine interaction of the order makes 
the baryon state masses mostly change by 1-25 MeV. The compositions and 
calculated masses are listed in Tables HI O [6l [71 [8], [9] and [TOl 

3. Results 

3.1. The ground-state baryons 

We can find that the calculated masses of the ground states agree well 



with the experimental values [35|. For the ground states, in addition, we 
compare our results of the order with the results using the order as in 
Table [21 It shows that the results of order are agree better with the 
experimental values, since most AM are somewhat smaller and the value of 



X;(AM)2 = 33 MeV^ is much smaller than that of 193 MeV^ in Ref. [31 

We further study the corrections to the Gell-Mann-Okubo mass formula 
(GMO) and Gell-Mann's equal spacing rule (GME). The corrections to Gell- 
Mann-Okubo mass formula for the baryon-octet and Gell-Mann's equal spac- 
ing rule for the baryon-decuplet can be written as 

^ = ^ + 5gmo, (40) 

Ms* -Ma = Ms* - Ms, + 5gmei 

= - Mb. + 5gme2, (41) 

Here as Sqmo = ^gmei = '^gme2 = 0, equations (140|) and pil) are the standard 
mass formulas which do not consider the -ffhyp term. However, in the real 
world Sqmo = —6.75 MeV, (5gmei = 8 MeV and Sqme2 = 11 MeV, which 
stand for the deviations of the above two mass formulas from experimental 
data. The deviations can be explained by a quark-mass dependent hyperfine 
interaction, together with the wave function size and mixing with excited 
states. In Tab. [3], we calculate and compare the deviations of the mass 
formulas. In the first row, it is the deviations obtained from experimental 
data. The deviations in the second row are results considering the one- 
gluon exchange process. After considering the higher order exchange process. 
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we list our deviations in the last row. It shows that the correction to the 
Gell-Mann-Okubo mass formula is distinctly improved and the correction to 
the Gell-Mann's equal spacing rule is also getting slightly better. Since the 
calculated H* mass is not very good, there is only little improvement on the 
corrections to the GME. However, the octet ground baryons are more close 
to reality and the corrections to the GMO improve significantly. It implies 
that the correction to the GMO is mainly due to the hyperfine interaction of 
Oial 



S ) 



3.2. Baryon excitations 

The excited baryon masses with spin J in the range 1/2-7/2 are listed in 
Tables m O El [HI El and [TUl We see that most of the masses of positive par- 
ity states are consistent with the experimental data. We compare the nucleon 
resonances iVi + (1440) and iVi-(1535). In particular the former is called the 

2 2 

Roper resonance and it is interpreted as the first positive parity excited state 
{iN'^S's^'^) or |A^2(56',0+)i"^)) with the iV = 2 harmonic oscillator band in 
the quark models, whereas the latter is interpreted as the first negative parity 
excited state ( \N^PmI~) or \N'^{70, with = 1. In the relativized 
quark model with one-gluon-exchange [36] or instanton-induced interactions 



37(1, the calculated mass of the iVi + (1440) exceeds the mass of the Ni- (1535) 

2 2 

by 80 MeV. It is not consistent with the experimental value that should be 
55-125 MeV below. In this work, however, using the nonrelativistic quark 
model with two-gluon-exchange, the calculated mass of A^i + (1440) is 1462 

2 

MeV, which locates below that of A^i - (1535), i.e., 1497 MeV. That can be ex- 

2 

plained as follows. First, in the harmonic-oscillator model (when U = and 
-f^hyp = 0), the nonstrange states of SU(6) multiplets with the same harmonic 
oscillator band are degenerate with E = {N + 3)a;. It makes the energy 
of the N = 2 positive parity excited states lu larger than that of the = 1 
negative parity excited states. However, the anharmonicity U = X]t<j Uijij) 
is a diagonal perturbation and break the initial degeneracy within the N = 2 
harmonic oscillator band. Thus the five SU(6) multiplets at = 2 have 
different energies, which we may denote by E(56',0"'"), £(70,0"*"), E(56,2+), 
E(70,2+) and E(20, 1+). The energy of (70,1") muhiplet at A^ = 1 is ex- 
pressed as E(70, 1~). These energies can be expressed by three parameters, 
i.e., a, b and Eq in ( I25i) -( l33l) . Using the reasonable parameters, E(56',0"'") 



becomes 1600 MeV, which is less than E(70, 1-)=1610 MeV [30|, ISlJ. That 



is to say that the baryon energies associated with the (56', 0"*") multiplet (as 
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indicated by the Roper resonance) are less than the baryon energies associ- 
ated with the (70, 1~) multiplet. Second, the hyperfine interaction further 
causes sphtting between the two states. It is noted that the diagonal matrix 
element of the contact term of 0{as) in the first positive parity excited state 
(iV2(56',0+)|^ I i/contact(0(a,)) | N\56',0+)1^) is -169 MeV, while that 
in the first negative parity excited state (A^^(70, 1^)| | -f^contact(0(Q;s)) | 
A^^(70, 1~)| ) is —67 MeV. The splitting between the two matrix elements 
is about —100 MeV, which further makes the energy of the first positive 
parity excited state lower. The contact term of 0{a^) makes the splitting 
between the two diagonal matrix elements increase by —7 MeV. In addition, 
the hyperfine interaction also brings the off-diagonal matrix elements, which 
not only cause mixing between the states with same A^, but also mixing be- 
tween the = and N = 2 band states. After diagonalizing the complete 
Hamiltonian matrices, we obtain the the baryon masses. The mass of the 
first positive parity state corresponding to A^i + (1440) is 1462 MeV. It is in 

2 

the experimental range of 1420-1470 MeV and is lower than that of the first 
negative parity corresponding to A^i-(1535). 

2 

Meanwhile, we find that the masses of the negative parity states are 
similar to the results of 0{as) [29]. Some predicated masses are consistent 
with the experimental data. However, some controversial states, such as 

Ai-(1405), A3- (1520) and A^i-(1535), are not consistent with the experi- 

22 2 

mental data. On one hand, the hyperfine interaction of O(a^) only makes 
the masses obtained from 0{as) change by 1-25 MeV. This effect is somewhat 
small and difficult to make these masses agree with the experimental data. 
On the other hand, these controversial states may have interpretations by in- 



troducing new physical contents beyond the 3q-quark model [38|, l39|, |40|, |41 



We find that the hyperfine interaction with two-gluon exchange does im- 
prove the calculated masses of some excited baryons, but not all of them. 
Thus it is still difficult to draw a conclusion that the quark model with the 
hyperfine interaction of O(a^) describes the excited baryon spectrum better 
than that of O(a^). 



4. Conclusions 

In this paper, we studied the baryon spectrum with hyperfine interaction 
of 0(ag) due to two-gluon exchange. We find that the matrix elements of 
hyperfine interaction with higher order are in the range (0.0014—0.0334)5, i.e. 
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the values of the matrix elements are in the range 0.4-10 MeV, which are one 
order less than those matrix elements 2-170 MeV obtained by leading order 
0(as). After diagonalizing the complete Hamiltonian matrices, the order 
makes the baryon state masses mostly change by 1-25 MeV. In addition, 
the order 0{a1) produces the calculated masses of ground states shift 7-10 
MeV comparing with that of the order 60 to 200 MeV generated by 0{as) in 
term of 5 = 270 MeV. The calculated masses are in good agreement with the 
experimental masses. The value ^(AM)^ shows that the predicted masses 
are closer to the experimental data by taking into account the higher order 
hyperfine interaction. In addition, the correction to the Gell-Mann-Okubo 
mass formula for the baryon-octet is distinctly improved and the correction 
to the Gell-Mann's equal spacing rule for the baryon-decuplet is also shghtly 
improved. 

It is noted that as becomes 0.7 when considering higher order hyperfine 



interactions, and it is smaller than the value 7r/2 in Ref. [3l[. This implies 
that the applicability of perturbative QCD calculations becomes more rea- 
sonable in the quark-model framework by taking into account higher order 
effects. In addition, we employ the same value of the parameter 6 for both 
ground and excited baryons. We assume that the value is about 270 MeV 
to consider the effects of the hyperfine interactions with higher order and 
the mixing between ground states and excited states. It is different from 
that the 6 value of 300 MeV for the excited states in Ref. Mi and 260 MeV 
for the ground states in Ref. 3ll|. Some parameters used in this paper are 



the same as those used in Ref. [30, |31|, 132|. Other parameters are roughly 



derived in terms of dynamical theory. These parameters can change in some 
reasonable ranges, but cannot change the order of magnitude of the physical 
quantities. 

Finally we conclude that the effects from higher order hyperfine inter- 
action should be considered in the quark model. Taking into account the 
hyperfine interaction of O(a^) due to two-gluon exchange, the calculated 
spectrum agree better with the experimental data. The higher order hyper- 
fine interaction bring also more realistic correction to the Gell-Mann-Okubo 
mass formula. Corrections due to two-gluon exchange on baryon excitations 
(including negative parity baryons) are also briefly discussed. 
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A. The hyperfine matrix elements in the SU(6) basis 

In this appendix we present detailed calculation of the hyperfine matrix 
elements of A in the SU(6) and uds bases. 

The spectrum of A, A and Q can be easily calculated in the SU(6) basis. 
However, E, A and S are usually first calculated by the matrix in the uds 
basis, which arc simpler. Then using the relation between the SU(6) and uds 
bases, we transform the matrix in the uds basis to the matrix in the SU(6) 
basis and obtain the baryon spectrum and the composition of the eigenstates 
in terms of the SU(6) basis. 

Take A for example, the SU(6) wave functions are 
(octet) 

= —-j={2uds — 2dus + usd — dsu — sud + sdu), 



(Al) 

$^ = --{usd- dsu + sud- sdu), (A2) 



(singlet) 



= — ^ [uds — dus — usd + dsu + sud — sdu) . 



The uds type wave function is 



and the spin wave functions are 



(A3) 



$A ^ -^{ud- du)s, (A4) 
V2 



X'+ = ^(tit-m), (A5) 
xl = ^(TU-iU), (A6) 



12 



X- 

xf/2 



-^(m + m -2 m), 



1 



(Tii + iti -2 iit), 



TTT, etc. 



(A7) 

(A8) 
(A9) 



The Mcis-type A states 







(AlO) 






(All) 


lA^-^pA^^) 




(A12) 



We display here only the top state of a given multiplet, and we also display 
only the state of maximum J from coupling a given L and 5"; the states of 
smaller J are constructed using standard tables in the LS order. 
The relation between the SU(6) and uds bases is that 



A^L 
A^L 
A^L 



AA 

PP 
p\ 



1 
V2 



w 

1 

2 



v 



^ ^ / 



Ai Lm 

As I'M 



(A13) 



where we use A instead of |A) to denote the state for simphcity. 



|A'L,,) 



IAk^L 



for L = 5" or D and 



lA^P, 



pX] 



lA^PpA) 



8 J-iM 



IAs'Pa), 
-lAi^P^). 



(A14) 



(A15) 
(A16) 



so the i^hyp matrix elements of the SU(6) and uds bases have the relation 



{asu{io)\H^^^\(3su{&)) = {oiuds\H^yl\(^uds) , 



uds I 



where 



trSU(6) 
-"hyp 



(A17) 

(A18) 
(A19) 
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Take the A with .P = 1/2+ with 
the SU(6) and uds bases is 



2 for example, the relation between 



/ A^S 



1 + 

AA2 



where 



A Spx^ 

A'/^pA^ 
A^PpAf 
V A^^pAf J 



5+ 





A 2q/ 1 + 


\ 




As 






Ai^^A/l 






As'^-Dm| 






A/P4^ 




v 


As^p^r 


/ 





1 


1 




2 


2 




1 


1 




2 


2 





+* 
































V 

The matrix i^contact in the -urfs basis with order a 



IS 



H. 



contact 



v 



_ 1 

2 



32 







X + 



X 





_3 
4 

3v5 
32 







(A20) 












\ 
























M 













-1 













-1 


/ 



(A21) 



7v^„ 
32 

32 











\ 













gX 




















































/ 



(A22) 



Using the relation in Eq. (]A2ip . the matrix H'^ 



contact 



contact 



contact 



in the SU(6) bases is 

(A23) 



/ 


_5 


-^(1 + a;) 


+ 4(1 -x) 











\ 






5 

16 


+ ^(l-x) 

















+^(l-a;) 


-Ul + 2x) 
























gX 
































\ 




















/ 
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In addition, we show some detailed calculation of the H^^ matrix elements. 
As mentioned in Eq. fl38l) . under exchange of quarks two and three, i.e., the 
relative coordinates become 

(23)p = Ip+ ^A, (A24) 
(23)A = ^p-^A, (A25) 



the spin functions are 



{23)x'=X, (A26) 

(23)x^ = + (A27) 

(23)x^ = ^x" - \x'- (A28) 
In term of the approximation ~ oa, the space functions are 

(23)Voo = ^00, (A29) 

(23)^00 = + (A30) 

(23)<o^ = \ro^ + \i^^^ + ^<. (A31) 

(23X = ^V,^-^i^l^ + \Vol (A32) 

(23)^^,^ = (A33) 

00 



where (23)V'^2^, (23)V'22 and (23)V'22 are similar with those of i2?>)^^^, (23)?/'^ 
and (23)^,"^ 



00 • 

Therefore, if |a) = = |A^5'pa1/2^), the matrix element -^contact is 

{xXv,o\HiL.a\xX<) = mr,i\x6\n,m?>)r,i) 

1 fiTrrv — * — * 
x^((23)xtl^i-^2|(23)x^) 

+ ^(^oo|5'(ri2)|^o'o)]-^ 

X = ^6. (A34) 



12v^m2 16 
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B. The hyperfine matrix elements with order ol 



The hyperfine matrix elements with order ctg have been calculated [29|, |30 
We calculate the matrix elements with order a1 by the same method. The 



simplest case corresponds to the interaction between quarks 1 and 2 [36|, |26 
The contact matrix elements are: 

>^{SaSm\Sl • S2\Sj3Sm)- 



9m? 



a. 



, 21a,;^2^ 1n(mV2p) +7j; ^,^ 



fBll 



The tensor term: 



HlLor = V{V2p)R2{l2) ■ S2(12), (B2) 



where R2(12) and 82(12) are the spatial-tensor and spin-tensor operators 
respectively. 

The calculation of the tensor matrix elements are done with the aid of the 



Wigner-Eckart theorem |42l . |43|, which is applied twice: once to the scalar 
product of the spin-tensor and spatial-tensor operators, and the other to the 
L = 2 tensor operator where L is made up from coupling Ip and l\. 

We apply the Wigner-Eckart theorem to the tensor product -f^tonsor 
obtain 

{a\Hll^J^) = {-iy''^--'m{L^LpS^Sp-2J) 



v/2L, + lv/25, + l(S„||S2(12)||5^) (B3) 
{L^npJpn^Jp^\V{^/2p)n2{l2)\\Lpnp^ 



The theorem is applied again to the spatial reduced matrix element 
{L^UpJp^nxJp^ I |1/(y2p)R2(12) \\Lpnp^ 
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xWilpJp^L^Lp; 2hJ^2L(, + 1^21,^ + 1 
x(n,J,J|l^(v^p)R2(12)|| 



where 



{rir, 



/,J|\/(v^p)R2(12)||n,,/,,) 
2(2L 



C(/,,200;/p,0) 



"PfS 

2L. + 1 



{'npJpAV{V2p)\np^^lp^;), 



and the hyperfine tensor energy V{-\/2p) with order a] is 



V{V2p) 



1 |4a, ^ 



3m2 (v^p)3 ' Stt 67r 
4 
3 



4 3q; 4 

x[ln(pGiiV^p) +7i^ - 7t] ^[In(mv^p) +7^ - -]}. 



TT 



(B4) 



(B5) 



(B6) 



For comparing with the hyperfine matrix elements with order a^, we also 
express the hyperfine matrix elements with order a1 in the unit of 5, which 
is 

An n"^ 

6 = ^ ^ ~ 270 MeV. (B7) 
3v27rm^ 



The 6 value of 300 MeV in Ref. [3(| is determined by considering only exited 
states with N = 2. When considering the effects of mixing between ground 



states with = and exited states with = 2 in Ref. 3l|], the 6 value 
is assumed to be 260 MeV. In this paper, we assume that this value differs 
slightly, which is 5 ~ 270 MeV, to consider the results of order effects and 
mixing between the ground states with = and the excited states with 
N = 2. 

In the calculation of the hyperfine matrix elements with order a^, there 
are more parameters needed. They are as, Pgr, ^d, a, nf and x, which 
satisfy 6 ~ 270 MeV. Since it is hard to derive the precise values of the 
parameters from foundational principle, we assume the reasonable values in 
term of the known ranges of the parameters. The light baryons belong to 
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the soft regime. It has been shown 4^, through a solution of Schwinger- 
Dyson equation, that when decreasing Q^, does not increase continuously 
but instead it saturates at a critical value 0^(0) ~ 0.7 ± 0.3. In our study 
regime, we assume as = 0.7. The effective mass range of the constituent 
quark is ~ 200 — 350 MeV, and we suppose rrid = 0.275 GeV here and 
the parameter a ~ 0.38 to satisfy the value of Eq. flB7p . By the similar 
relation that fXcGR oc rric for charmonium in Ref. [26[, we suppose that the 
renormalization scale fidGR = 0.333 GeV for baryons with mi = m2 = md, 
and /isGR = 0.333/x GeV for baryons with mi = m2 = m^, where we take 
X = 0.616. The number of effective quark flavors is Uf = 2, which is the 
number of quarks with mass less than the energy scale /igr- 

In summary, we use the values = 0.7, x = 0.616, a = 0.38, = 
0.275 GeV, fidCR = 0.333 GeV, /x.gr = 0.333/a; GeV and 6 = 270 MeV. The 
hyperfine matrix elements are all listed in Tables [11] and [T2l 
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Table 1: The confinement energy in MeV. 



Confinement energy 





5 = 


S = -l 


S = -2 


S = - 


Emj 

E[^L] 


1135 


1296 


1457 


1618 


1610 


1770 


1825 


1985 


1610 


1700 


1895 


1985 


Em^] 


1705 


1895 


1910 


2070 


E[roo] 


1810 


1945 


2040 


2150 


E^] 


1705 


1805 


2000 


2070 




-105 


-100 


-90 


-85 


E[<J 


2020 


2145 


2225 


2315 


EW.I] 


1890 


2085 


2055 


2215 


E[r2t] 


1935 


2065 


2150 


2245 




1890 


1975 


2175 


2215 




-40 


-40 


-35 


-35 
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Table 2: The ground-state baryons (in units of MeV). 



State (J 


^) Mass 


Mass [31] 




AM 


AM [31] 




939 


940 


939 





-1 




1113 


1110 


1116 


3 


6 




1192 


1190 


1193 


1 


3 


^1 + 


1317 


1325 


1318 


-1 


-7 


Af 


1231 


1240 


1232 


1 


-8 


V* 3 + 
^ 2 


1383 


1390 


1385 


2 


-5 


■^*3 + 


1526 


1530 


1530 


4 







1673 


1675 


1672 


-1 


-3 




where the value AM 


= Mass 


-Mexp 
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Table 3: The corrections to GMO and GME in MeV. 





^GMO 


^GMEl 


'^GME2 


^exp 


-6.75 


8 


11 


m 


2.5 


10 


5 




-4.75 


9 


5 
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Table 4: The calculated N spectrum and composition in the SU(6) basis. 



State 


Mass 


Exp. 




Composition 










(MeV) 


(MeV) 














Nl+ 

2 


939 


939 


0.94 


0.25 


0.22 


-0.01 


0.00 " 




2 


1462 


1420 - 1470 


-0.26 


0.96 


0.04 


0.02 


0.00 




2 


1748 


1680 - 1740 


0.19 


0.08 


-0.90 


0.37 


0.09 




]Vi + 


1887 


1885 - 2125 


-0.06 


-0.05 


0.37 


0.83 


0.40 




2 
2 


2060 




0.01 


0.01 


-0.07 


-0.41 


0.91 


N^PAt 


2 


1734 


1700 - 1750 


" -0.16 


0.86 


-0.49 


0.01 


0.00 " 




2 


1857 


1855 - 1975 


-0.78 


-0.30 


-0.25 


0.46 


-0.13 




iV3 + 
2 


1952 




0.57 


-0.07 


-0.29 


0.64 


-0.41 




2 


1975 




0.16 


-0.41 


-0.77 


-0.03 


0.31 




2 


2058 




0.10 


0.07 


0.10 


0.51 


0.84 




2 


1738 


1680 - 1690 


0.88 


-0.47 


0.00 ■ 








2 


1953 




-0.45 


-0.84 


-0.31 








2 


2007 




-0.15 


-0.27 


0.95 








2 


1943 


1940 - 2090 1.00 













2 


1497 


1525 - 


1545 


0.80 0.59 " 


N^Pm\ 


2 


1650 


1645- 


1670 


-0.59 0.80 _ 


N^Pmf 


2 


1548 


1515- 


1525 


■ 0.99 -0.13 ■ 


N^Pm\- 




1731 


1650 - 


1750 


0.13 0.99 


N'PmI 




1655 


1670- 


1680 


1.00 


N^Pm\ 
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Table 5: The calculated A spectrum and composition in the SU(6) basis. 



State Mass Exp. Composition 



(MeV) (MeV) 





1867 






" 0.56 


0.83 " 










1900 


1870- 


1920 


0.83 


-0.56 










1231 


1232 


" 0.98 


-0.21 


0.00 


0.00 " 




A|+ 


1790 


1550- 


1700 


0.21 


0.96 


0.13 


-0.08 




Af+ 


1904 


1900- 


1970 


0.03 


0.15 


-0.92 


0.37 




Af 


1972 






0.00 


0.02 


0.37 


0.93 




Af 


1931 


1865- 


1915 


0.92 


0.39 " 








Af 


1967 


1724- 


2200 


-0.39 


0.92 










1902 


1915- 


1950 1.00 












1668 


1600- 


1660 


1.00 








A'Pm\~ 


Af 


1668 


1670- 


1750 


1.00 











Table 6: The calculated Q spectrum and composition in the SU(6) basis. 



State 


Mass 
(MeV) 


Exp. 
(MeV) 






Composition 








2182 






" 0.61 


0.79 " 










2202 






0.79 - 


-0.61 










1673 


1672 




0.99 


-0.15 0.00 


0.00 " 






2078 






-0.15 


-0.98 -0.08 


0.04 




"2 


2208 






0.01 


0.08 -0.96 


0.27 




"2 


2263 






0.00 


0.02 0.28 


0.96 






2224 






0.97 


0.23 " 










2260 






-0.23 


0.97 








at 


2205 




1.00 










a\- 


2015 






1.00 












2015 






1.00 
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Table 7: The calculated A spectrum of positive parity. 



State 


Mass 


Exp. 




Composition 














(MeV) 


(MeV) 


















Ai+ 

■"■2 

Ai+ 


1113 


1116 


0.95 


-0.24 


-0.17 


0.05 


0.01 


0.00 


0.00 " 




1606 


1560 - 1700 


-0.22 


-0.95 


0.09 


-0.19 


-0.02 


0.00 


0.00 


Ai+ 

■"■2 


1764 


1750 - 1850 


0.12 


0.18 


0.15 


-0.96 


0.02 


0.00 


0.00 


A 20 1 + 


Ai+ 


1880 




0.16 


0.02 


0.91 


0.16 


-0.34 


-0.02 


-0.08 




Ai+ 


2013 




0.05 


-0.01 


0.33 


0.08 


0.85 


0.08 


0.39 


A *n 1 + 


Ai+ 


2173 




0.01 


0.00 


0.06 


0.02 


0.40 


-0.28 


-0.87 


A8^P^i+ 


Ai+ 


2198 




0.00 


0.00 


0.00 


0.00 


0.04 


0.96 


-0.29 


A|+ 

A3 + 


1836 


1850 - 1910 


0.11 


-0.80 


-0.34 


0.47 


-0.01 


0.07 


0.00 ■ 


A *<? 3 + 


1958 




0.58 


-0.04 


0.63 


0.31 


-0.21 


-0.34 


0.06 


Af 


1993 




-0.70 


-0.28 


0.48 


0.09 


0.24 


-0.36 


-0.06 


Al^/^Mf 


A3 + 


2061 




0.38 


-0.10 


0.00 


-0.21 


0.80 


-0.12 


-0.39 


As^Dsf 


A|+ 


2121 




-0.12 


0.49 


-0.11 


0.77 


0.19 


0.07 


-0.30 


A *n 3+ 


A3 + 


2134 




-0.02 


-0.14 


0.48 


-0.02 


0.00 


0.83 


-0.24 


A/P^f 


■"■2 

Af 


2174 




_ -0.04 


-0.07 


-0.09 


-0.16 


-0.47 


-0.20 


-0.83 


As^P^f 


A5 + 
■"■2 


1839 


1815 - 1825 


" -0.82 


-0.33 


0.47 


0.00 


-0.02 " 








A|+ 


2008 




-0.16 


0.91 


0.35 


0.06 


0.13 








A5 + 


2103 




-0.43 


0.18 


-0.61 


-0.62 


0.14 






As^Z^mI* 


A|+ 


2129 




-0.35 


0.12 


-0.52 


0.73 


-0.24 






As^Z^Aff 


Af 


2155 




0.02 


0.12 


0.08 


-0.27 


-0.95 










2064 


2020 - 2140 1.00 
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Table 8: The calculated S spectrum of positive parity. 



State 



K 

^2 
V3 + 
^2 
V3 + 

E- 

V3 + 

^2 
V3 + 
^2 

^2 

V5 + 

? + 
^2 



Mass 
(MeV) 



Exp. 
(MeV) 



Composition 



V 20; 1 + 

V 2c. 1 + 

^10 5m _ 

Ti 4 rr 3 + 

V 4 a/ 3 + 

V 4 Q 3 + 

V 4r) 3 + 



1192 


1193 


0.98 


-0.14 


-0.16 


0.02 


0.00 


0.01 


0.00 ■ 






1664 


1630 - 1690 


-0.17 


-0.96 


-0.19 


0.08 


0.00 


0.00 


0.00 






1924 


1826 - 1985 


-0.12 


0.17 


-0.86 


-0.20 


-0.18 


0.36 


0.07 






1986 




0.04 


-0.08 


0.35 


0.03 


-0.81 


0.44 


0.12 






2022 




-0.02 


0.10 


-0.09 


0.94 


0.09 


0.28 


0.10 






2069 




-0.03 


0.08 


-0.24 


0.26 


-0.54 


-0.67 


-0..36 






2172 




-0.01 


0.02 


-0.06 


0.02 


-0.10 


-0.38 


0.92 






1383 


1385 


0.98 


0.17 


0.01 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


1868 


1720 - 1925 


-0.16 


0.92 


-0.33 


0.02 


0.07 


-0.04 


0.05 


-0.04 


0.00 


1947 




-0.01 


0.07 


0.09 


-0.94 


-0.01 


0.22 


0.07 


0.20 


-0.01 


1993 




0.03 


-0.13 


-0.46 


-0.05 


-0.69 


-0.28 


0.44 


0.12 


-0.08 


2039 




0.05 


-0.26 


-0.62 


-0.20 


0.54 


-0.43 


-0.17 


-0.04 


0.04 


2075 


2070 - 2140 


-0.03 


0.13 


0.53 


-0.13 


0.05 


-0.79 


0.20 


-0.12 


-0.07 


2098 




-0.01 


0.08 


0.06 


0.06 


-0.21 


-0.25 


-0.55 


0.71 


0.27 


2122 




0.00 


0.02 


-0.03 


-0.20 


-0.39 


-0.08 


-0.46 


-0.64 


0.42 


2168 




0.00 


-0.02 


0.03 


0.06 


0.18 


0.05 


0.47 


0.09 


0.86 


1949 


1900 - 1935 


" -0.95 


0.04 


0.20 


-0.03 


0.22 ' 










2028 




-0.04 


-0.88 


0.13 


0.44 


-0.08 










2062 


2051 - 2070 


-0.21 


-0.15 


-0.97 


-0.03 


0.02 










2107 




0.21 


-0.04 


-0.02 


0.11 


0.97 










2154 




-0.05 


0.44 


-0.09 


0.89 


-0.08 











1+ 

2 



Mr, 



D 

^8 J^S2 



^2 


2002 


2025 - 2040 


■ 0.87 


-0.50 " 




^9. 


2106 




0.50 


0.87 J 
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Table 9: The calculated S spectrum of positive parity. 



State 



-1 + 

=!+ 

'2 
-1 + 
"2 
-1 + 
■"1 
-1 + 

% 

"2 
-3 + 
'2 
-3 + 
'2 
-3 + 

= 1+ 

= 1+ 

= 1+ 

= 1+ 
'2 
-3 + 
'2 
-6 + 
'2 
-6 + 



-5 + 
'2 
-5 + 
■"1 



Mass Exp. 
(MeV) (MeV) 



Composition 



:^ 2c 1 + 

"8 '-'SJ 

2 c/ 1 + 
"8 ■5^2 
r-. 2c 1 + 

" 2 c 1 + 
"10 »5M2 

= 10 -t^S^ 

- 4n 1 + 

"8 -f^Mj 

- 2 p 1 + 
"8 t^A2 

- 4 c 3 + 
"10 S'S^ 

4 c' 3 + 

"10 

^ 4c 3 + 
=8 '-'M2 

2 n 3 + 

=8 Us 2 
= 4 n 3 + 

"10 

~ 2n 3 + 

=8 

= 4n 3+ 

=8 

" '^n 3+ 

=10 

2p 3 + 
"8 i^A2 

=8 ■L'Sa 



1317 
1750 
1982 

2054 
2107 
2149 
2254 

1526 
1952 
1970 
2065 
2114 
2174 
2184 
2218 
2252 

1959 
2102 
2170 
2205 
2239 



1318 



1530 



-0.96 


0.22 


0.17 


0.02 


0.00 


-0.01 


0.00 ■ 






0.21 


0.97 


-0.08 


0.09 


0.00 


0.03 


0.00 






-0.17 


-0.05 


0.88 


0.03 


0.18 


0.38 


0.08 






-0.06 


0.13 


-0.38 


-0.11 


-0.72 


-0.54 


-0.16 






-0.01 


-0.08 


-0.05 


0.97 


0.05 


-0.21 


-0.09 






0.02 


-0.03 


0.16 


0.20 


-0.65 


0.59 


0.38 






-0.01 


0.01 


-0.06 


-0.01 


0.14 


-0.41 


0.90 






-0.98 


-0.17 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


-0.07 


0.40 


0.28 


-0.71 


0.02 


0.47 


-0.03 


-0.14 


0.00 


0.14 


-0.82 


-0.27 


-0.42 


-0.07 


0.23 


0.01 


-0.04 


0.00 


0.03 


-0.14 


0.44 


0.15 


-0.63 


0.13 


-0.56 


0.15 


0.11 


0.05 


-0.33 


0.74 


0.16 


0.54 


0.08 


0.12 


0.00 


-0.03 


0.01 


-0.05 


0.06 


-0.09 


0.05 


-0.40 


-0.31 


-0.83 


0.19 


-0.01 


0.05 


-0.32 


0.25 


0.46 


0.45 


-0.52 


0.04 


0.38 


0.00 


0.01 


-0.07 


0.42 


-0.18 


0.54 


0.17 


-0.49 


-0.47 


0.00 


-0.01 


0.07 


0.14 


-0.24 


0.18 


0.53 


-0.14 


0.76 



0.84 
0.03 

-0.21 
0.46 

-0.19 



0.00 
0.86 
0.14 
-0.18 
-0.46 



-0.52 0.00 0.15 

0.08 0.49 0.12 

-0.57 0.04 -0.78 

0.58 0.32 -0.57 

-0.25 0.81 0.19 



= 4n 5 + 

"10 Um2 



- 2 n ."i + 

"8 J^M2 

= 4n + 

"8 Um2 

= 2 n 5 + 

"10 Um2 



= 7 + 
"2 


2074 


■ 0.81 


0.58 ' 


4n 7+ 

"10 Us 2 


= 7 + 
"2 


2189 


0.58 


-0.81 


=10 Um2 
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Table 10: The calculated spectrum of the negative-parity S=-l and S=-2 excited baryons. 



State 


Mass 


Exp. 




Composition 






(MeV) 


(MeV) 










Af 


1559° 


1402 - 1410 


" 0.09 


0.50 


0.86 ■ 






1656 


1660 - 1680 


0.59 


0.66 


-0.45 




H' 


1791 


1720 - 1850 


0.79 


-0.56 


0.24 






1560" 


1520 


" -0.02 


0.47 


0.88 






1702 


1685 - 1695 


0.13 


-0.87 


0.47 


A 2p 3- 




1859 




0.99 


0.12 


-0.05 


A 2p 3- 




1803 


1810 - 1830 


1.00 








^' 


1657 


1620 


0.69 


-0.71 


-0.15 




^' 
^' 


1746 


1730 - 1800 


0.70 


0.60 


0.38 




1802 




-0.18 


-0.37 


0.91 




^2 


1698 


1165 - 1685 


" -0.15 


-0.95 


-0.29 


n V ''p 3- 




1790 




0.89 


-0.01 


-0.45 


V 2 p 3 — 


^' 


1802 




0.42 


-0.33 


0.85 


V 2 p 3 — 


% 


1743 


1770 - 1780 


1.00 








"2 


1772 




" 0.48 


0.85 


-0.19 ' 


4p 1~ 


"2 


1894 




0.75 


-0.29 


0.59 


2p 1- 
"8 -fjlf2 


"2 


1926 




0.44 


-0.43 


-0.78 


^ 2 p 1- 

=10 -nv/2 


= 3- 
"2 


1801 


1817- 1828 


-0.1 


0.95 


-0.28 " 


"Ip 3~ 
=8 -0/5 


= 3- 
"2 


1918 




-0.09 


0.27 


0.96 


2p 3- 
=8 -fM2 


"2 


1976 




0.99 


0.12 


0.06 


= 2p 3- 
"10 ^M2 


"2 


1917 




1.00 






= 4p 5- 
—8 "Af2 



" Tliese results do not consider the second-order effects in Ref. [29] 
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Table 11: Matrix elements i?contact(0(Q;s)) (in units 



CO 
o 



/ 0.0293 -0.014 

-0.0014 0.0057 

-0.0020 -0.0334 

0.0000 0.0000 

\ 0.0000 0.0000 

/ 



0.0000 0.0000 \ 
0.0000 0.0000 



-0.0020 
-0.0334 

0.0183 

0.0000 

0.0000 0.0000 0.0309 / 



0.0000 0.0000 
-0.0259 0.0000 



A 10 S'g 



3 + 



/ -0.0293 0.0020 

0.0020 -0.0057 

0.0000 0.0000 

0.0000 0.0000 



0.0000 0.0000 

0.0000 0.0000 

-0.0208 0.0000 

0.0000 -0.0259 



0.0293 


0.0014 


0.0010 - 0.0059X 


-0.0010 - 0.0059X 


0.0000 


0.0000 


0.0000 \ 


0.0014 


0.0057 


-0.0167 - 0.0147X 


0.0167 - 0.0147X 


0.0000 


0.0000 


0.0000 


0.0010 - 0.00591 


-0.0167 - 0.0147X 


0.0028 + 0.0206X 


-0.0028 + O.OI5O1 


0.0000 


0.0000 


0.0000 


-0.0010 - 0.0059.T 


0.0167 - 0.01471 


-0.0028 -1- O.OI5O1 


0.00284 + 0.0094X 


0.0000 


0,01)1)0 


0,0001) 


0.0000 


0.0000 


0.0000 


0.0000 


-0.0103 - 0.01841 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


-0.0103 - 0.0206X 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


-0.0103 + 0.0413X J 


-0.0098 -1- 0.0541X 


0.0005 - O.Ollli 


0.0003 - 0.00981 


-0.0003 - 0.0020a: 


0.0000 


0.0000 


0.0000 \ 




0.0005 - O.Ollli 


-0.0019 + 0.02631 


-0.0056 - 0.0245a: 


0.0056 - 0.0049X 


0.0000 


0.0000 


0.0000 




0.0003 - 0.0098.T 


-0.0056 - O.O24.5.T 


0.0145 + 0.01321 


0.0164 -0.0188.T 


0.0000 


0,0000 


0.0000 




-0.0003 - 0.0020.T 


0.005<i - 0.1)049.1 


0.0164 - O.OI881 


0.0145 + (l.0244,E 


0.0000 


0,01)00 


0.0000 




0.0000 


0.0000 


0.0000 


0.0000 


-0.0069 - 0.01761 


0.0028 - 0.00471 


0.0000 




0.0000 


0.0000 


0.0000 


0.0000 


0.0028 - 0.0047X 


-0.0069 - 0.0199X 


0.0000 




0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





0.0309 I 







-0 


0098 - 0.0270i 


-0.0005 + 0.00663! 


-0.0005 - 0.00281 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 \ 


-0 


0005 + 0.00561 


-0.0019 - 0.0132a: 


0.0079 - 0.0069X 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 




-0 


0005 - 0.0028X 


0.0079 - 0.0069X 


-0.0019 -O.OI881 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 






0.0000 


0.0000 


0.0000 


-0.0070 + 0.0353X 


0.0000 


-0.0020 - 0.0166a: 


0.0000 


0.0020 - 0.0033X 


0.0000 






0.0000 


0.0000 


0.0000 


0.0000 


-0.0070 - 0.0176X 


0.0000 


0.0028 - 0.0047X 


0.0000 


0.0000 






0.0000 


0.0000 


0.0000 


-0.0020 - 0.0166a: 


0.0000 


0.0120 + 0.0176a: 


0.0000 


0.0189 - 0.0199X 


0.0000 






0.0000 


0.0000 


0.0000 


0.0000 


0.0028 - 0.0047X 


0.0000 


-0.0070 - 0.0199X 


0.0000 


0.0000 






0.0000 


0.0000 


0.0000 


0.0020 - 0.0033X 


0.0000 


0.0189 - 0.0199X 


0.0000 


0.0120 + O.O22I1 


0.0000 


V 




0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0309 / 



Table 12: Matrix elements i?tensor(0(Q;^)) (in units 



CO 



Ai'fUr 

V 4q 3 + 
4 0/3 + 

^10 -i^A/o 



/ 0.0000 
0.0000 
0.0000 
0.0000 

\ 0.0000 

/ 0.0000 
0.0000 
0.0000 
0.(100(1 
0.0000 
0.0000 

\ 0.0000 

/ 0.0000 
0.0000 
0.0000 
0.1)00(1 
0.0000 
0.0000 

\ 0.0000 

/ 0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 

\ 0.0000 



0.0000 
0.0180 
0.0008 
0.0187 
-0.0204 



0.0000 0.0000 

0.0000 0.0000 

0.0000 0.0000 

0.0180 0.0008 

0.0000 0.0000 

0.0000 
0.0000 

0.0000 
0.00(11) 

0.01531 0.0090 - 0.0037^ 



0.0000 
0.0000 

0.0000 
(1.0000 



0.0000 
0.0000 



0.0000 
0.0000 



0.0000 
0.0000 
0.0000 
0.0(100 
0.0120 - 0.01021 
0.0120 + O.OOSlx 
0.0000 



0.0000 N 
0.0000 
0.0000 
-0.0204 
0.0000 / 

0.0000 
0.0000 

0.0000 
0.00(11) 
0.0090 - 0.01431 
0.0000 
0.0000 

0.0000 
0.0000 
0.0000 
0.0000 

-0.0085 + 0.0072a; 0.0085 - 
-0.0085 + 0.0178X 0.0085 ■ 
0.0000 



Aio S'g 



S2 



'St , 



/ 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0255 -0.0180 

0.0000 0.0255 0.0000 0.0205 

\ 0.0000 -0.0180 0.0205 0.0000 



0.0000 
0.01531 

0.0090 - 0.00371 
0.0090 - 0.0143.C 
-0.0100 - 0.02701 



0.0000 
0.0000 

0.0000 
0.0000 
-0.0096 + 0.01513: 



-0.0096 -I- 0.01511 0.0072 -I- 0.0226i 
-0.0136 - 0.0107X -0.0101 + O.OI6O1 



0.0000 \ 
0.0000 

0.0000 

o.oooo 

-0.0130 - 0.01071 

-0.0101 + o.oiaox 
0.0000 



0.0000 0.0000 0.0000 0.0000 \ 

0.0000 0.0120-0.01021 0.0120 + 0.00511 0.0000 

0.0000 -0.0085 + 0.00721 -0.0085 + 0.0178i 0.0000 

0.0000 (1.0085 - 0.01441 0.0085 - 0.0072a: 0.0000 

0.01441 -0.0097- 0.01681 -0.0097 + 0.01131 0.0000 

0.0072a: -0.0097 + O.OllSx -0.0097 - 0.0352x -0.0321x 

0.0000 0.0000 -0.0321a: 0.0000 J 



0.0000 

0.0000 

0.0000 
-0.0086 + 0.0072X 
-0.0170 + 0.01441 
0.0060 - 0.0051a: 
-0.0170 - 0.0072a; 
-0.0060 - 0.0102X 

0.0000 



0.0000 

0.0000 

0.0000 
-0.0085 - 0.0036X 
-0.0170 - 0.0072» 
0.0060 - 0.0126X 
-0.0170 + 0.0250a: 
-0.0060 + 0.0051X 

0.0000 



0.0000 
-0.0086 + 0.0072X 
-0.0085 - 0.0036X 

0.0000 
0.0097 - 0.0113X 

0.0000 
0.0097 + 0.0057a: 

0.0000 

0.0000 



0.0000 
-0.0170 + 0.0144X 
-0.0170 - 0.0072X 
0.0097 - 0.0113X 

0.0000 
-0.0068 + 0.0080X 

0.0000 
0.0068 + 0.0160X 

0.0000 



0.0000 
0.0060 - 0.0051X 
0.0060 - 0.0126X 

0.0000 
-0.0068 + 0.0080X 

0.0000 
-0.0068 + 0.0127X 
0.0000 
0.0000 



0.0000 
-0.0170 - 0.0072X 
-0.0170 + 0.0250X 
0.0097 + 0.0057X 

0.0000 
-0.0068 + 0.0127X 

0.0000 
0.0068 - 0.0080X 

0.0000 



0.0000 0.0000 \ 

-0.0060 - 0.0102X 0.0000 
-0.0060 + 0.0051X 0.00000 



0.0000 
0.0068 + 0.0160X 

0.0000 
0.0068 - 0.0080X 

0.0000 

0.0000 



0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 / 



